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ABSTRACT 

The possibility of determining MDM model parameters on the basis of observable 
data on the Abell-ACO power spectrum and mass function is analysed. It is shown 
that spectrum area corresponding to these data is sensitive enough to such MDM model 
parameters as neutrino mass m u , number species of massive neutrino N u , baryon content 
Qb and Hubble constant h = Hq/ 100km/ s/Mpc. The x 2 minimization method was used 
for their determination. If all these parameters are under searching then observable 
data on the Abell-ACO power spectrum and mass function prefer models which have 
parameters in the range Q u (~ 0.4 — 0.5), low Qf, (< 0.01) and h (~ 0.4 — 0.6). The 
best-fit parameters are as follows: N v = 3, m v = 4AeV, h = 0.56, Qt, < 0.01. The 
high-fit ~ 0.4 — 0.5 solutions are obtained when mass of neutrino is fixed and < 3eV. 

To explain the observable excessive power at k ~ 0.05h/Mpc the peak of Gaussian 
form was introduced in primordial power spectrum. Its parameters (amplitude, position 
and width) were determined along with the MDM model parameters. It decreases x 2 > 
increases the bulk motions, but does not change essentially the best-fit MDM parameters. 

It is shown also that models with the median Q u ~ 0.2 — 0.3 {m v ~ 2.5, N v ~ 2 — 3) 
and Qf, = 0.024/ h 2 , which match constraints arising from cosmological nucleosynthesis 
and high redshift objects, are not ruled out by these data (A% 2 < 1). 

Subject headings: Large Scale Structure: Abell-ACO power spectrum, mass function, 
Mixed Dark Matter models, initial power spectra, best-fit cosmological parameters 
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1. Introduction 

The observable data on large scale structure of the 
Universe obtained during last years and coming from 
current experiments and observational program give 
a possibility to determine more exactly the param- 
eters of cosmological models and the nature of the 



dank matter. Up till now the most certain data are 
about the largest sCale inhomog'eneities of the Cur- 
rent particle horizon of the order of ~ 7000/i _1 Mpc 
(h = f/o/100 km/s/Mpc, Hq is today Hubble con- 
stant) which are obtained from the study of all-sky 
temperature fluctuations of cosmic microwave back- 
ground (CMB) with ~ 10° angular resolution by the 



space experiment COBE (3moot ct al. 1992, Ben 
ictt ct al. 1994, Bennett et al. 1996). According to 



them the primordial power spectrum of density fluctu- 
ations is approximately scale invariant P p i(k) = Ak n 
with n — 1.1 ± 0.2 that well agrees with the predic- 
tions of standard inflation model of the Early Uni- 
verse (n — 1, Slo = !)■ Besides, they most certainly 
determine the amplitude of a linear power spectrum 
(or normalization constant A) which does not depend 
on any transition processes, nonlinearity effects and 
other phenomena connected with the last stages of 
large scale structure formation. On the contrary, the 
CMB temperature fluctuations at degree and sub- 
degree scales as well as the space distributions of the 
cluster of galaxies, galaxies, quasars, Lyman-a clouds, 
etc. are defined by those processes and also depend 



ess entially on the nature of the dark matter. — Theo- 



reti cally it is taken into account by introducing the 

transfer function T(k) which transforms the primor- 
dial (post-inflation spectrum) into the postrecombi- 
nation (initial) one - P(k) = P pi (k)T 2 (k), which de- 
fines all characteristics of the large scale structure of 
the Universe. The transfer function depends also on 
the curvature of the Universe or the present energy 
density in units of critical density, Qq, vacuum en- 
ergy density or cosmological constant content of 
baryons fif,, and values of the Hubble constant. The 
theory of a large scale structure formation is so far 
advanced today that all these dependencies can be 



acq uiaLelv calculated foi the fixed model by public 



ava ilable codes (e.g. — CMBfast one by Beljak fc Zal- 



iar riaga 199G )r — The actual problem now is the dc 



As most advanced candidates for the dark matter 
are cold dark matter (CDM), particles like axions, 
hot dark matter (HDM), particles like massive neu- 
trinos with m v ~ 1 — 20eV and baryon low luminosity 
compact objects. The last ones can not dominate as 
it results from the cosmological nucleosynthesis con- 
straints { Q b h 2 < 0.024, [Tytler et al. 1996| , |Songalia et 



al. 1997, Schramm and Turner 1997) and observation 



of microlensing events in the experiments like MA- 
CHO, DUO, etc. The pure HDM model conflicts with 
the existence of high redshift objects, the pure CDM 
one, on the contrary, overpredicts them. Therefore 
mixed dark matter model (CDM+HDM+baryons) 
with Qhdm = &v < 0.3 looks more viable. The 
advantage of these models is a small number of free 
parameters. But today it is understood already that 
models with the minimal number of free parameters, 
such as a standard cold dark matter (sCDM, one pa- 
rameter) or a standard cold plus hot mixed dark mat- 
ter (sMDM, two parameters) only marginally match 
the observable data. A better agreement between the- 
oretical predictions and observable data is achieved in 
the models with a larger number of free parameters 
(tilted CDM, open CDM, CDM or MDM with the cos- 



mological term, see review in Valdarnini et al. 1998 
and references therein). 

The oscillations of solar and atmospheric neutri- 
nos registered by SuperKamiokande experiment show 
that the difference of rest masses between r— and 
/i-ncutrino s is 0.02 < Am T fl < 0.08eU |Fukuda et 



al. 1998, Primack & Gross 199?. It also gives a 
lower limit for the mass of neutrino m v > |Am| and 
does not exclude models with cosmologically signif- 
icant values ~ 1 — 20eU. Therefore, at least two 
species of neutrinos can have approximately equal 
masses in this range. Some versions of elementary 
particle theories predict m„ e m Ur ~ 2.5eV and 
m v „ w m u . ~ 10~ 



denote the electron 
accordingly (e.g 



eV , where v e , v T , and v s 



and sterile n eutrinos 
Berezhiani et al. 1995|) . The 



strongest upper limit for the neutrino mass comes 
from the data on a large sca le structure of our Uni- 
y, m„j93h 2 < 3 (|Holtzman 1989|, |Davis 



ter mination of the nature of the dark matter and the 

rest of the above mentioned parameters by means of 
comparison of theoretically predicted and observable 
characteristics of the large scale structure of the Uni- 



Summers fc Schlegel 1992, Schaefer fc Shafi 1992, Van 



Dalen fc Schaefer 1992|, Novosyadlyj 1994, |Pogosyan 
fc Starobinsky 1995| , pVla. 1996| , |Valdarnini et al. 1998| ) 



that for h — 0.8 (the upper observable limit for h) 
gives y i m n < 18eV. It is interesting that the up- 
per limit for the mass of electron neutrino obtained 
from supernova star burst SN1987A neutrino signal 
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is approximately the same m Ve < 20eV. 

Is it possible to find the best fit neutrino mass 
from experimental data on a large scale structure of 
the Universe? The problem is that it must be de- 
termined together with other large number uncer- 



tain parameters such as h, £l m , Qb, etc. Here we 



study the possibility of finding them by \ mini- 
mization method. Realization of such a task became 
possible in principle after the appearance in litera- 
ture of accurate analytical approximations of trans- 
fer function for mixed dark matter model in at least 
4-dimension space of the above mentioned cosmolog- 



ical parameters T(k; f^, m V) N u , h) (Eiscnstcin & Hu 
19SfrR [Novosyadlyj et al. 199§|). That is why that 



even CMBfast codes are too bulky and slow yet for 
sea rching the cosmological parameters by the meth- 



ods of minimization of x , like Levenberg-Marquardt 



one (see Press et al. 1992 ). 

The next problem is a choice of the observable data 
suitable for the solution of this task. They must be 
enough accurate, sensitive to those parameters and 
not too dependent on the model assumptions about 
the formation and nature of objects. The most sensi- 
tive to the presence of neutrino component are scales 
of order and smaller of its free-streaming (or Jeans) 
scale k > kj(z) = 8 (j^f) /VTT^h^ 1 Mpc be- 
cause perturbations at these scales are suppressed and 
it is imprinted in the transfer function of the HDM 
component. At z ~ for cosmologically significant 
neutrino masses it is approximately galaxy clusters 
scale. The power spectrum reconstructed from space 
distributions of galaxies is distorted significantly by 
nonlincarity effects the accounting of which is model 
dependent (Peacock & Dodds 1994). The models of 



the formation of smaller scale structures or high red- 
shift objects (e.g. Lyman-a damped systems, Lyman- 
a clouds, quasars etc.) contain the additional as- 
sumptions and parameters which makes their using 
rather problematic in such an approach. The CMB 
temperature anisotropy at subdegree angular scales 
(first and second acoustic peaks) has minimal addi- 
tional assumptions (e.g. secondary ionization) but its 
sensibility to the presence of neutr ino component is 
low (< 10%, Dodelson et al. 1995 ). These data are 



sensitive and suitable for determination by x mini- 
mization methods other set of parameters such as tilt 
of primordial spectrum n, f2o, h, f^, £l\ or/and pa- 
rameters of scaling seed models of structure formation 



see Lineweaver fc Barbosa 1997, Durrer et al. 1997) 



tion mass of rich clusters seem to be suitable for de- 
termining the best fit values of ra„ and N v because 
they do not depend on above mentioned additional 
assumptions. 



The data on rich clusters power spectrum ( Einasto 



et al. 1997) were used by Eisenstein et al. 1997 
and Atrio-Barandela et al. 1997 for analyzing ~ 
100/i -1 Mpc clustering. The first collaboration group 
tried to explain the narrow peak in the power spec- 
trum at ~ 100ft. -1 Mpc scale by baryonic acous- 
tic oscillations in low- and high-Oo models (Oq = 
&CDM + ^b)- In both cases such an approach needs 
very high content of baryons fib (> 0.3), that is es- 
sentially out of the cosmological nucleosynthesis con- 
straints. The second one has shown that this feature 



is in a greement with Saskatoon data (|Ncttcrficld et al 



1997 ) on AT/T power spectrum at subdegree angu- 



lar scales. They have concluded that these data pre- 
fer models with built-in scale in the primordial power 
spectrum which can be generated in the more com- 
plicated inflation scenario (e.g. double one). 

For reducing the number of free parameters we re- 
strict ourselves to analysis within the framework of 
the matter dominated Universe and standard infla- 
tion scenario: Q m = Qq = 1, n = 1 without the 
tensor mode of cosmological perturbations. The free 
parameters in our task will be baryon content Of,, di- 
mensionless Hubble constant h, neutrino mass m v , 
and number species of neutrinos with equal masses 
N v . 

The outline of this paper is as follows: the observ- 
able data which will be used here are described in 
Section 2. The method of determination of parame- 
ters and its testing are described in Sect. 3. Results 
of best fit finding of parameters under different com- 
bination of free and fixed ones are presented in Sect. 
4. Discussion of results and conclusions are given in 
Sect. 5 and 6 accordingly. 

2. Experimental data set 

The most favorable data for the search of best fit 
cosmological parameters are real power spectrum re- 
constructed from redshift-s pace distribution of Abell 
ACQ clusters o f galaxies ( |Einasto et al. 1997 , Ret 



The data on Abell- ACO power spectrum and func- 



zlaff et al. 1997). It is biased linear spectrum reliably 
estimated for 0.03 < k < 0.2h/Mpc whose position 
of maximum (k max w .05h / M pc) , inclination before 
and after it are sensitive to baryon content Ob, Hub- 
ble constant h, neutrino mass m„ and number species 
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of massive neutrinos N v (see Fig. 1-4). Here in numer- 
ical calculations the data of last estimation of power 
spectrum by Retzlaff ct al. 1997 will be used. All 



Table 1: Experimental data set. 



the sources of systematic and statistical uncertainties 
as well as window function and differences between 
Abell and ACO parts of sample have been accurately 
taken into account there. The values of the Abell- 
ACO power spectrum for 13 values of k PA+ACo(kj) 
(j = 1, 13) and their la errors are presented in Table 
1 and are shown in figures. 

Other observable data which will be used here 
are constraints of amplitude of the fluctuation power 
spectrum at cluster scale derived from cluster mass 
and X-ray temperature functions. It is usually formu- 
lated as a constraint for density fluctuations in top- 
hat sphere of 8/i _1 Mpc radius, erg, which can be easy 
calculated for the given initial power spectrum P{k): 



a l = ^> \ fc2p ( fc ; n b, h, m v , N v )W 2 (8k/h)dk, 



1 

(1) 

where W(x) = 3(sinx — xcosx)/x 3 is Fourier trans- 
formation of top-hat window function. The different 
collaboration groups gave similar results which are in 
the range of a s ~ 0.5 — 0.7. The new optical de- 
termination of the mass function of nearby galaxy 
clusters (|Girardi et al. 1998| ) gives median values: 
a § = 0.60 ± 0.04[j. It match es very well the cluste r 
X-ray temperature function ( Viana fc Liddlc 1996 ). 



For taking into account the data of other authors I 
shall be more conservative and will use it with 3tr er- 
ror bars instead of la one. But, as we will see, it does 
not rule out predicted a$ value from the la limit of 
the observable one by Girardi et al. 199S for best fit 
parameters determined here. 

The COBE 4-year data will be used here for nor- 
malization of power spectra. A useful fit for them is 
the amplitude of density perturbation of the horizon 
crossing scale Sh, which for a flat model with the n = 1 
eq uals S h = 1.94 • 10" 5 ( |Liddle et al. 1996j , [Bunn and 
Wlite 1997). Taking into account the definition of 



<5/j (Liddle et al. 1996) and the power spectrum, the 
normalization constant A is calculated as 

A = 277 2 Sl(3000/h) 4 Mpc 4 . 



1 Jeans scale for neutrino component in all cases analysed here 
is smaller than the cluster scales therefore all the matter is 
clustered and the term f2 0,6 in the original form is omited 



IV", . 

INO 


k 


Vi 




Ay, 


i 


0.030 


9.312 


10 4 


±59723.65 


2 


0.035 


1.037 


10 5 


±65488.2 


3 


0.040 


1.039 


10 5 


±58014.15 


4 


0.047 


1.258 


10 5 


±51005.75 


5 


0.054 


1.448 


10 5 


±68638.6 


6 


0.062 


1.016 


10 5 


±39184.6 


7 


0.072 


8.098 


10 4 


±25179.7 


8 


0.083 


5.444 


10 4 


±21925.45 





0.096 


5.303 


10 4 


±24914.75 


10 


0.11 


3.853 


10 4 


±13344.5 


11 


0.13 


2.031 


10 4 


±8546.35 


12 


0.15 


2.039 


10 4 


±9804.3 


13 


0.17 


1.691 


10 4 


±9383.21 


14 




0.60 




±0.12 



3. Method and its testing 

The Abell- ACO power spectrum is connected with 
matter one by means of the cluster biasing parameter 
b c r- 

PA+Aco{k) = b 2 cl P{k;fl b ,h,m u , N u ). (2) 

For fixed parameters fib, h, m v , N v and b c i the val- 
ues of PA+ACo(kj) are calculated for the same kj as 
in Table 1 and as according to (1). Let's denote 
them by y 3 (j = 1, ...,14), where j/i , yi 3 corre- 
spond PA+Aco(ki), ■-,PA+Aco(ki3), and y u is a 8 . 
Their deviation from observable data set (noted by 
the tilde) can be described by x 2 : 



14 

E 



Vi ~ Vi 



(3) 



where jjj and Ajjj are experimental data set and their 
dispersion accordingly. Then parameters Qb, h, m„, 
N v and b c i or some part from them can be determined 
by minimizing % 2 using Levenberg-Marquard method 
(Press et al. 1992). The derivatives of predicted val- 
ues on search parameters which are required by this 
method will be calculated numerically. The step for 
their calculation was experimentally assorted and is 
10~ 5 of the values for all parameters. 

The analytical approximation of MDM transfer 
function will be used in the form: 

T M DM(k; fib, h, m v , N u ; z) = 

T C DM+b(k; fib, h; z)D(k; fl b , h, fl u , N v ] z), 
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where TcpM+b{k\ 0;,, h; z) is the transfer function 
by Eisenstcin & Hu 1997a for CDM+baryon sys- 



tem (z is redshift), the correction factor for the 
HDM component D(k) was used in the form given 
by Novosyadlyj et al. 1998. It is correct in a suffi- 



ciently wide range of search parameters (for a more 



et al 



detailed analysis of its accuracy see in Novosyadlyj 
1998j ) . We suppose the scale invariant primor- 



dial power spectrum because the initial power spec- 
tra of MDM models now is as follows: Pmdm (&) = 
AkTlj DM (k; Q b , h, m v ,N v ;z). 

The method was tested in the following way. I 
calculated the MDM power spectrum for the given 
parameters (e.g. fi& = 0.15, fl u = 0.2, N u = 1, 
h = 0.5) using CMBfast code, normalized to 4-year 
COBE data, calculated as and interpolated P(k) for 
the same kj (j — 1, 13) which are in Table 1. Then 
I have took cluster biasing parameter b c i = 3 and cal- 
culated model PA+Aco(kj)- The 'experimental' er- 
rors for them as well as for as I have suggested to be 
the same as relative errors from Table 1. These model 
experimental data like the ones in Table 1 were used 
for search of parameters fif,, h, f2„, and b c i (N u is 
fixed and the same). The initial (or start) values of 
the parameters I have put as random deviated from 
the given ones. In all cases the code found all the 
given parameters with high accuracy. 

4. Dependence of density fluctuations power 
spectra at cluster scale on cosmological pa- 
rameters 

Before finding of the best-fit parameters let's look 
how the power spectrum of density fluctuations at 
cluster scale depends on search parameters. For this 
we leave only b c i as a free parameter and fix the re- 
maining ones. In Fig.l such a dependence of rich 
cluster power spectra on is shown for h = 0.5, 
fib = 0.05 and N u = 1. The r.m.s. of density fluc- 
tuations in the top-hat sphere of 8/i _1 Mpc radius 
in models with f2„ =0.1, 0.2, 0.3, 0.4 are a 8 = 0.93, 
0.81, 0.75, 0.71 accordingly. The best-fit values of b c i 
are presented in the caption of Fig.l. The deviations 
of the predicted rich cluster power spectra and mass 
function in these models from the observable ones are 
correspondingly x 2 = 17.3, 9.88, 6.64, 5.33. There- 
fore, for the MDM model with h = 0.5, Q b = 0.05 
and N v = 1 Abell-ACO power spectrum and mass 
function prefer high f2„ (~ 0.3 — 0.4). 

Now we repeat the same calculations for different 




k [h Mpc"'] 

Fig. 1. — The rich cluster power spectrum for MDM 
models with different f2„ (N u , Qb and h are fixed). 
The filed circles are experimental Abell-ACO power 
spectrum by Retzlaff et al. 1997. The best-fit biasing 
parameters for the models with fl u =0.1, 0.2, 0.3, 0.4 
are b c i =2.6, 2.8, 2.9, 2.9 accordingly. 

number species of massive neutrinos N v = 1, 2, 3 
and fixed £1„ = 0.2 (Fig. 2). The as's for these 3 mod- 
els are 0.81, 0.73, 0.68 accordingly, the corresponding 
deviations of predicted rich cluster power spectra and 
mass functions from the observable ones respectively 
are x 2 =9.88, 6.48, 5.54. So, the MDM model with 
three species of equal mass neutrino is preferable. 

In the first two cases (h fixed and equal 0.5) the 
mass of neutrino was different for differing Sl u (N v 
fixed) and N v {VL V fixed) because they are connected 
by relations 

m v = mci v h 2 /N u . (4) 

Let's fix the neutrino mass [m v = 2.5eV), suggest 
that N v — 2 and repeat calculations for different 
h =0.5, 0.6, 0.7. The results are shown in Fig. 3. <r$ 
for these 3 models are following 0.71, 0.98, 1.24. The 
X 2 for all points of power spectrum and as are 5.72, 
19.9 and 42.6 accordingly. Therefore, when neutrino 
mass is fixed (by laboratory experiments for example) 
the data prefer low h. 

Similarly, one shall calculate rich cluster power 
spectra for different f2& when the rest of the param- 
eters are fixed. The results for f2{, =0.05, 0.1, 0.15, 
0.2, 0.25, 0.3 are presented in Fig. 4. The correspond- 
ing cr 8 's are following 0.71, 0.64, 0.58, 0.53, 0.48, 0.44, 
the characteristics of deviations of the predicted val- 



5 




Fig. 2. — The rich cluster power spectrum for MDM 
models with a varying number of species of massive 
neutrino N v (Vt v , fib and h are fixed) . The filed circles 
are the same as in Fig.l. The best-fit biasing param- 
eters for models with N v = 1, 2, 3 are b c i =2.8, 3.1, 
3.3 accordingly. 

ues from the observable ones x 2 f° r these models are 
5.72, 4.28, 3.61, 3.70, 4.59, 6.39. The minimum X 2 is 
for model with SI5 = 0.15. 

As we see the theoretically predicted values of the 
chosen data are sensitive to search parameters m„, 
AT„, and h. It is interesting now where the global 
minimum of \ 2 m space of these parameters is when 
all or a part of them are free. 

5. Results 

The searching oim u , N v , fib and h by \ 2 Levenberg- 
Marquardt minimization method can be realized in 
the following way. We shall put m„, Ot, , h and b c i 
or part of them free and find the minimum of \ 2 for 
N v =l, 2, 3 in a series. The lowest value from them 
will be suggested as minimum of x 2 for each set of 
free parameters. This is because the N v possesses the 
discrete value. 

The key point is narrowing the range of search pa- 
rameter values. The analytical approximation of the 
MDM power spectra used here is accurate enough in 
the following range of parameters: 0.3 < h < 0.7, 



n„ < 0.5, tt b < 0.3, N v < 3 flNovosyadlyj et al. 1998|) . 
By the upper and lower boundaries of h, fi„ and fib 
availability of the used analytical approximation we 



Fig. 3. — The rich cluster power spectrum for MDM 
models with a varying h (m v , N v and fib are fixed). 
The filed circles are the same as in Fig.l. The best- fit 
biasing parameters for models with h =0.5, 0.6 and 
0.7 are b c i =3.2, 2.5, 2.1 accordingly. 



admeasure the range of search values of these param- 
eters. We make these boundaries as 'mirror walls'. 

5.1. All parameters are free 

The minima of \ 2 m a 4-dimensional space of pa- 
rameters fi„, f2b, h and b c i for models with 1, 2 and 3 
species of massive neutrinos are achieved for the set 
of parameters presented in Table 2. The spectra for 
them are shown in Fig. 5 and cts's are presented in the 
Table 2. (The accuracy of analytical approximation 
of MDM spectra is better than 5%). 

As we can see \ 2 is f ew times lower than the formal 
degree of freedom, d = n — m, where n is the number 
of data points, m is the number of free parameters. 
The reason is that not all the points of the Abell- 
ACO power spectrum presented in Table 1 are inde- 
pendent. The numerical experiment has shown that 
the minimal number of points which determine the 
same MDM parameters is w 7 (odd points of Pa^) 
in Table 1, for example). Indeed, such a spectrum can 
be described by amplitude and inclination at small 
and large scale ranges and the second order curve at 
the peak (or maximum) range. It means that real 
d fa 3-4. 

Therefore, in the 5-dimension space of free parame- 
ters (f2„, N v , fib, h and b c i) the global minimum of \ 2 
is achieved for the MDM model with 3 sorts of mas- 
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Fig. 4. — The rich cluster power spectrum for MDM 
models with a varying Ot, (r?v, N v and h are fixed). 
The filed circles are the same as in Fig.l. The best- 
fit biasing parameters for models with il b =0.05, 0.1, 
0.15, 0.2, 0.25 and 0.3 are b cl =3.2, 3.4, 3.7, 3.9, 4.1 
and 4.3 accordingly. 

Table 2: Best-fit parameters of MDM models with 1, 2 
and 3 sorts of massive neutrinos for Abell-ACO power 
spectrum by Retzlaff et al. 1997 and mass function 
by Girardi et al. 1998. 
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sive neutrinos. It has the lowest m„ and the highest h 
which better matches the data on immediate measure- 
ments of Hubble constant. However, it is unexpected 
that the found Cl u is so high and fl b is so low. They 
contradict the data on high redshift objects and nucle- 



osy nthesis constraint (0.007 < fl b h 2 < 0.024, [Tytler 
etfl. 1996, 3ongalia et al. 1997, Schramm and Turner 



1992)" accordingly. The MDM models with so high a 
Q„ (~ 0.4 — 0.5) also have a problem with the galaxy 
formation, <7q ~ 1 for them. Let's analyze the cases 
with additional constraints which can lead us out of 
this difficulty. 



Fig. 5. — The rich cluster power spectra of MDM 
models with best-fit parameters f2„, Sit, , h and b c i 
for 1, 2 and 3 sorts of massive neutrinos (Table 2). 
The filed circles are experimental Abell-ACO power 
spectrum by Retzlaff et al. 1997. 

5.2. Coordination with nucleosynthesis con- 
straint 

The increasing of baryon content can decrease this 
difficulty (see Eisenstein et al. 1997). We shall fix 
baryon content by the upper limit which is resulted 
from the nucleosynthesis constraint Q, b h 2 = 0.024 and 
keep up the rest parameters as free. The found best- 
fit parameters are in the Table 3, rich power spectrum 
for the case with 3 sorts of massive neutrino is shown 
in Fig. 6 (dotted line). The spectra for the cases with 
1 and 2 sorts are close to this one. As we can see fl u 

Table 3: Best-fit parameters of MDM models with 1, 2 
and 3 sorts of massive neutrinos for Abell-ACO power 
spectrum by Retzlaff et al. 1997 and mass function 
by Girardi et al. 1998 when baryon content is fixed 
by nucleosynthesis constraint (il b h 2 = 0.024). 
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increases when Of, decreases and the minima of x 2 are 
achieved at high again. But they are quite close 
to the corresponding minima from the previous table 
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(Ax 2 < 1) 



,— , 'CT r 



' C ' r 




k [h Mpc~'] 

Fig. 6. — The rich cluster power spectrum of MDM 
models with 3 sorts of massive neutrinos and best- 
fit parameters for the cases when all parameters are 
free (solid line), when baryon content fib is fixed by 
nucleosynthesis constraint (dotted line), when mass 
of neutrino m u = 2.5eV is fixed (dashed line) and 
when both fib and m„ are fixed (dashed dotted line). 
The filed circles are the same as in Fig. 5. 

5.3. When the mass of neutrino is known 

An interesting question ensuing from last two items 
is: which best-fit values of fib and h can be obtained 
from these data on the Abell-ACO power spectrum 
and mass function in the case when mass of neutrino 
is determined by any physical or astrophysical experi- 
ments and is known. Let's assume that m„ is fixed but 
the number of species N v is unknown. We fix fl v by 
relation (4) and the rest of parameters leave free. The 
search in such an approach was unsuccessful because 
it halted in the upper limit of f2b=0.3. When this 
'mirror wall' was removed the solutions were found 
but with extremely high content of baryons for which 
an accuracy of analytical approximation for MDM 
spectra is worse (~ 15 — 20%). Results for m v —2.beV 
and 3eV are presented in Table 4. The rich cluster 
power spectrum for m v = 2.5eV and N v = 3 is shown 
in Fig. 6 (dashed line). The spectra for 1 and 2 sorts 
are close to this one. 

The x 2 's in all cases here are lower than in Table 2 
because the performance of analytical approximation 
of MDM spectra for so high a fib and h is essentially 



Table 4: Best-fit parameters of MDM models with 1, 2 
and 3 sorts of massive neutrinos for Abell-ACO power 
spectrum by Retzlaff et al. 1997 and mass function 
by Girardi et al. 1998 when neutrino mass is fixed 
(m„=2.5 and 3.0eV, fl v = m v N v /92,h 2 ). 



N v 


x 2 


fl v (m„) 


fi b 


h 


hi 


^8 


1 


1.53 


0.05 (2.5) 


0.47 


0.75 


3.22 


0.65 


2 


1.44 


0.09 (2.5) 


0.45 


0.79 


3.21 


0.65 


3 


1.39 


0.12 (2.5) 


0.43 


0.82 


3.20 


0.65 


1 


1.51 


0.06 (3.0) 


0.47 


0.75 


3.21 


0.65 


2 


1.42 


0.10 (3.0) 


0.44 


0.79 


3.20 


0.65 


3 


1.37 


0.14 (3.0) 


0.42 


0.83 


3.19 


0.65 



worse than in the allowance range. Therefore we can 
not conclude that the global minimum of x 2 hi the 
4-dimension space of parameters m v , N v , fib and h is 
in the range of high fib and h. It is in point with the 
parameters which are in the last row of Table 2. But 
we certainly conclude that when m u ~ 2 — 3eV the 
minimum is absent in the range of fib < 0.3, 0.3 < h < 
0.7. Therefore the Abell-ACO power spectrum and 
mass function among the MDM models with m„ < 
AeV and N v < 3 prefer fl b > 0.3 and h ~ 0.8 that 



agrees well with the results by Eiscnstein et al. 1997 



5.4. fib and m v are fixed 

One can look now which h is preferable by Abell- 
ACO power spectrum and mass function when neu- 
trino mass and baryon content are fixed by the other 
observable constraints or theoretical arguments. Let's 
put that m v = 2.5eV (f2„ = m v N v /9?>h 2 ) and fib = 
0.024/ft 2 is fixed by the upper limit of nucleosynthesis 
constraint. Only h and b c i are free parameters. Their 
best- fit values found for 1, 2 and 3 sorts of massive 
neutrino are presented in the Table 5. The rich clus- 
ter power spectrum for N v — 3 MDM model with 
those parameters is shown in Fig. 6 (dashed dotted 
line). The spectra for 1 and 2 sorts are close to this 
one. 

As we see in the MDM model with 3 sorts of 2.5eV 
neutrinos the best-fit value of h and ct§ are closer 
to the corresponding observable data than in mod- 
els with 1 or 2 sorts. 
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Table 5: Best-fit parameters of MDM models with 1, 2 
and 3 sorts of massive neutrinos for Abell-ACO power 
spectrum by RetzlafF et al. 1997 and mass function by 
Girardi et al. 1998 when baryon content and neutrino 



mass are fixed: 
m u N v /93h 2 ). 



fl h = 0.024//i 2 



=2.5eV 



K 


x 2 


Q v (m„) 




h 


b c i 


^8 


1 


4.20 


0.16 (2.5) 


0.15 


0.41 


3.90 


0.55 


2 


3.31 


0.24 (2.5) 


0.11 


0.47 


3.76 


0.56 


3 


2.85 


0.29 (2.5) 


0.09 


0.52 


3.69 


0.56 



6. Discussion 

Rich cluster power spectra of models with the best 
fit parameters are within the error bars of the cor- 
responding experimental data (Fig. 5-6). But none of 
them explains the peak at k » 0.05h/Mpc that cor- 
responds to the linear scale sa 120/i _1 Mpc . It has 
excess power at ~ 50% in comparison with the best- 
fit model and ~ 30% in comparison with the high-f^ 
on e. It is mo re prominent yet in the data by Einasto 
et d. 1997 . Apparently, it is a real feature of the 
power spectrum. The necessity of a similar feature 



in the power spectrum was argued earlier by the ex- 
pla nation of Great Attracto r phenomenon ( Hnatyk et 
al. |l995j iNovosyadlyj 1996|) . A s ample of the Abcll- 
ACO clusters of galaxies used by RetzlafF et al. 1997 
is placed in 60° double-cone with the axis pointing to- 
wards the Milky Way pole. The Great Attractor, on 
the contrary, placed in the plane of our galaxy. There- 
fore, they are an independent experimental demon- 
stration of the reality of those peak. Other impor- 
tant arguments for its validity come from pencil-beam 
redshift survey by Broadhurst et al. 1990 and from 



2-dimensional power spectrum of the Las Campanas 
Redshift Survey (Landy et al. 1996). The angular 



correlations in the APM survey (Gaztanaga & Baudh 
19S0) and high-rcdshift absorption lines in quasar 
spectrum (Quashnock et al. 1996) also show similar- 



features at these scales. It was shown also by Atrio- 
Bai fandcla et al. 1997j that this ~ 120/i _1 Mpc peak 
well agree with Saskatoon data on the AT/T power 
spectrum. Therefore, the data used here on rich clus- 
ter power spectrum are based on the surveys which 
represent a fair sample of ~ 120/i -1 Mpc structures 
and that peak is significant despite the large error 
bars of experimental data. 



Obviously, that turnabout to open (il < 1) mod- 
els or flat with cosmological term (fio + Q\ = 1) 
does not improve the situation with the explanation 
of that peak in our approach. It is because the max- 
imum of power spectra in those models is shifted to 
larger scales in comparison with matter dominated 
flat models analyzed here. Explaning of it by baryonic 
acoustic oscillations calls for extremely high content 
of baryons that disagree with nucleosynthesis con- 
straint (see Eisenstein et al. 1997). Therefore we face 
a necessity to consider models with a built-in scale in 
the primordial power spectrum again. 

Let's determine the parameters of this peak. The 
comparison of rich cluster power spectrum predicted 
by the MDM model with the best-fit parameters 
(Table 2) with the observable one showed that the 
peak has approximately the Gaussian form. There- 
fore we approximate it by the function p(k) = 1 + 
a p exp(2(k p — k) 2 /w 2 ), where a p , k p and w p are am- 
plitude, center and width of the peak accordingly. We 
set the power spectrum in the form of PMDM+p(k) = 
Pmdm (Jk\ Hfe, h, m v , N u )p(k; a p , k p , w v ), and repeat pre- 
vious calculations with additional free parameters a p , 
k p and w p . 

It should seem that this peak causes such high best- 
fit values of f2„ or in Tables 2-4. The results of the 
search for best-fit parameters in the 8-dimensional 
space of the MDM+peak model parameters showed 
that it is not so, that well agrees with the numerical 

results by Retzlaff et al. 1997. The introducing of the 

.2 



peak really decreases the \ but the MDM model pa- 
rameters are changed weakly. It is because they are 
determined mainly by the inclination of the Abell- 
ACO power spectrum after the peak and as as the 
most accurate value of the data set used here. The 
models with 3 sorts of massive neutrino are preferable 
like in the previous cases. In Table 9 the best-fit pa- 
rameters of the MDM models with 3 sorts of massive 
neutrino as well as best-fit parameters of the peak are 
presented for 4 cases: all the MDM parameters were 
free (1st row), baryon content Qb was fixed by the up- 
per limit of nucleosynthesis constraint (2), neutrino 
mass was fixed at m„ = 2.5eV (3), Qb and ra v were 
fixed (4). The x 2 for them are 0.81, 0.86, 1.11, 1.04 
accordingly. In all the cases except (3) the as = 0.6, 
in (3) case the as = 0.66. The rich cluster power 
spectrum for these cases are shown in Fig. 7. 

The introducing of such a peak increases the pre- 
dicted bulk velocities in a top-hat sphere of the radius 
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Table 6: Best-fit parameters of MDM+peak mod- 
els with 3 sorts of massive neutrinos for Abell-ACO 
power spectrum by RetzlafF et al. 1997 and mass func- 
tion by Girardi et al. 1998. The fixed parameters are 



noted by W (fi 6 = 0.024/h 2 



=2.5eV). 



rn v 


fib 


h 




kp 


a p 


uu p 


4.6 


0.01 


0.58 


3.14 


0.056 


0.46 


0.011 


5.0 


0.064^*) 


0.61 


3.11 


0.056 


0.47 


0.012 


2.5W 


0.424 


0.82 


3.16 


0.054 


0.34 


0.007 


2 rW 


084« 


n 53 


3 33 


060 


63 


013 



k [h Mpc~'] 

Fig. 7. — The rich cluster power spectrum of 
MDM+peak models with 3 sorts of massive neutrino 
and best-fit parameters from Table 6. The filed circles 
are the same as in Fig. 5. 



R whose r.m.s. values can be calculated according to 

1 

dkP MD M{k)W 2 (kR), 



V 2 - ¥1. 



where W(kR) is the Fourier transform of this sphere. 
So, for R = 50/i _1 Mpc it increases from 340fcm/s 
to 360fcm/s for the best-fit model (3rd row of Table 
2, 1st row of Table 6) and from 330fcm/s to 345km /s 
for a model with fixed m v and fif, (3rd row of Table 
5, the last row of Table 6). The observable value 
of bulk velocity for this scale is V50 = 375 ± 8 5fcm/s, 
wh ich follows from Mark III POTENT results ( |Kolatt 
& Dekel 1997). Therefore, this peak is preferable also 
by the data on large scale peculiar velocity of galaxies 
and Great Attractor like structures. However, the 
models with high values of Vt w ~ 0.4 — 0.5 (m v ~ 



4 - 7eV), which are best-fit ones for the Abell-ACO 
data, have problems with the explanation of galaxy 
scale structures and high redshift objects. But models 
with median Q„ ~ 0.2 - 0.3 (m u ~ 2.5, N v ~ 2 - 3) 
are not ruled out by these data (A\ 2 < 1). On the 
contrary, the CDM model with f2f, < 0.2 and h > 0.5 
is ruled out by these data at a high confidence level 
because for them A% 2 < 15. 

At last it must be noted that primordial spectrum 
feature like this peak is inherent for double inflation 



models (Kofman et al 
Kofman 



1985, Kofman & Linde 1987 



Pogosyan 1988, Gottloeber et al. 1991, 



Polarski & Starobinsky 1992) and inflationary model 




wherein an inflation field evolves through a kink in the 
potential ( Btarobinsky 199^ ) . Both classes of these 
models were confronted with the observational data 



on the A bell- ACO power spectru m by Lesgourgues et 
al. 1997| and |Rctzlaff et al. 1997| accordingly. 



Conclusions 



The Abell-ACO power s pectrum by [Retzlaff et aT 



1997 and mass function by Girardi et al. 1998 in the 
parameter space of the MDM model (Slo = 1) prefer 
a region with high fl u (~ 0.4 - 0.5), low Q b (< 0.01) 
and h (~ 0.4 — 0.6). The best-fit parameters are as 
follows: N v = 3, m v = 4.4eV, h = 0.56, fi 6 < 0.01. 
Unfortunately, experimental uncertainties of the data 
used here for the determination of these parameters 
give no chance to rule out models with a different set 
of parameters at a sufficiently high confidence level. 
The MDM models with baryon content at the upper 
limit of the nucleosynthesis constraint (fl^h 2 = 0.024) 
do not outstep A\ 2 = 1 of best-fit model (see Table 
3). The high-fib (~ 0.4 — 0.5) solutions are obtained 
when neutrino mass are fixed and < 3eV. 

Introducing artificially into the primordial power 
spectrum a peak of Gaussian form decreases the x 2 , 
increases the bulk motions but does not change es- 
sentially the best-fit parameters of the MDM models. 
It means that determinative for these parameters is 
mainly inclination of the Abell-ACO power spectrum 
at the scales smaller than the scale of the peak posi- 
tion and (jg as the most accurate value of the data set 
used here. 

Hereby, the power spectrum of the Abell-ACO 
clusters of galaxies and mass function are a sensitive 
test for the MDM model parameters. But more accu- 
rate data on power spectrum of matter density fluc- 
tuations are necessary for more certain determination 
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of cosmological parameters. 
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